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Abstract 

Recent observations of Type la supernovae provide evidence for the acceleration of our universe, 
which leads to the possibility that the universe is entering an inflationary epoch. We simulate 
it under a "big bounce" model, which contains a time variable cosmological "constant" that is 
derived from a higher dimension and manifests itself in 4D spacetime as dark energy. By properly 
choosing the two arbitrary functions contained in the model, we obtain a simple exact solution 
in which the evolution of the universe is divided into several stages. Before the big bounce, the 
universe contracts from a A-dominated vacuum, and after the bounce, the universe expands. In the 
early time after the bounce, the expansion of the universe is decelerating. In the late time after the 
bounce, dark energy (i.e., the variable cosmological "constant") overtakes dark matter and baryons, 
and the expansion enters an accelerating stage. When time tends to infinity, the contribution of 
dark energy tends to two third of the total energy density of the universe, qualitatively in agreement 
with observations. 

PACS numbers: 04.50. +h, 98.80.-k, 98.80.Es. 
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I. INTRODUCTION 



Recently, more and more observations of Type la supernovae suggest that the observable 
universe is presently undergoing an accelerating expansion, ^'^ which is contrary to what has 
always been assumed that the expansion is slowing down due to gravity with a positive 
decelerating parameter q. This, together with the recent high precision measurement of 
the cosmic microwave background (CMB) fluctuations,^'^ indicates the existence of dark 
energy characterized by a negative pressure, contributing with about 70%^'^'^'^ of the totel 
energy density of the universe (the others being essentially dark matter and baryons) . The 
simplest model for dark energy is the cosmological constant A. However, because of the 
so-called cosmological constant problem,''''^ one wishes to have a dynamically decaying A or 
an evolving large-scale scalar field known as quintessence^'^ to explain the acceleration of 
our universe. 

It is of great interest that our conventional universe is embedded in a higher-dimensional 
world as required in the Kaluza-Klein theories and the brane world theories. In this paper, 
we consider a five-dimensional cosmological model presented by Liu and Wesson. ^° Rather 
than the "big bang" singularity of the standard cosmology, this 5D model is characterized 
by a "big bounce", at which the "size" of the universe is finite. Before the bounce the 
universe contracts and after the bounce the universe expands. This model is 5D Ricci-flat, 
implying that it is empty viewed from 5D. However, as is known from the induced matter 
theory,^^'^^ 4D Einstein equations with matter could be recovered from 5D Kaluza-Klein 
equations in apparent vacuum. This approach is guaranteed by Campbell's theorem that 
any solution of the Einstein equations in N- dimensions can be locally embedded in a Ricci- 
fiat manifold of (N+1)- dimensions. An important result of the 5D bounce model is that 
a time variable cosmological "constant" can be isolated out, in a natural way, from the 
induced 4D energy-momentum tensor. 

The 5D bounce solution contains two arbitrary functions ^{t) and ^{t). It was shown in 
Ref. [10] that by properly choosing //(i) and v{t), one can obtain exact solutions suitable 
to describe both the radiation-dominated universe and the matter-dominated universe as in 
the standard FRW models. In this paper we will show that by properly choosing ii(t) and 
z/(t), we can also obtain exact solutions suitable to describe our present accelerating universe 
which is dominated by dark energy. 
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II. AN ACCELERATING UNIVERSE MODEL 



The 5D cosmological solution reads^° 

dS^ = B^de - { + r'dnA - dy\ 
y 1 — kr'^ ) 

fi ot n 

Here /i = yu(t) and z/ = v{t) are arbitrary functions, /c is the three-dimensional curvature 
index {k = ±1,0), and is a constant. This solution satisfies the 5D equations Rab = 
{A,B = 0123; 5). So one has R = and R^^Rab = 0. The 5D Kretschmann invariant is 
found to be 

I = RabcdR'^^'^'^ = • (2) 

which shows that K determines the curvature of the five-dimensional manifold. This solution 
was firstly derived by Liu and Mashhoon in a different notation in Ref. [14]. The bounce 
property was firstly discussed by Liu and Wesson. Further studies concerning the bounce 
singularity and other properties can be found in Ref. [15-20]. This solution can also be used 
to construct exact brane cosmological models. 

The 4D part of the empty 5D metric gives an exact solution of the 4D Einstein 
equations with an effective or induced energy-momentum tensor. It was shown in Ref. [10] 
that this energy-momentum tensor could be modeled by a perfect ffuid with density p and 
pressure p, plus a cosmological term A: 

^'^^T„/3 = (p + p)UaUf3 + (A - p)ga(3- (3) 

where u" = dx'^/ds is the 4-velocity. Furthermore, suppose the equation of state being of 
the form p = '-fp, then the solution gives 

2 ( p? + k pp\ 



1 + 7 V AA, 



1 + 7 



1 + 37 ( p^ + k\ pp 



(4) 
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We can see that if the two functions and z/(t) are given, then the two scale factors A(t, y) 
and B{t,y) are fixed by (Q) and the mass density p{t,y) and the cosmological term A{t,y) 
are also fixed by (jH). Generally speaking, on a given y = const hypersurface, A = A(t) is a 
time-variable cosmological term. 

From the metric ^ we see that the form Bdt = (^A/ fij dt is invariant under an arbitrary 
coordinate transformation t t{t). This would enable us to fix one of the two arbitrary 
functions fi{t) and leaving another to account for the various content of the cosmic 
matter. Since A is a variable function, one can not simply set A = or A = constant 
without losing the generality of the solution. Therefore, to compare the 5D solution with 
the A = FRW models, a constraint was used in Ref. [10] that the cosmological term A 
shown in (0]) should decay faster than the density p decreases at late-times of the universe. 
Under this constraint, it was found that the choice fi{t) oc t~^^^ and z/(t) oc t^^^ with k = 
and K = 1 gives a matter-dominated model, and the choice p{t) oc t^^/^ and u = const 
with k = and K = 1 gives a radiation-dominated model. For these two cases, the universe 
evolves, at late-times, with the same rate as in the FRW models. 

Here we are going to compare the 5D solution with the present accelerating universe. We 
are aware that neither one of the two terms A and p in Q is negligible presently. Therefore 
we require that A should decay not faster than p does in the late-times. In what follows we 
will show how this is achieved. 

Astrophysical data are compatible with a. k = flat 3D space. Meanwhile, we assume 
the density p being composed mainly of cold dark matter and baryons. So we choose A; = 0, 
K = 1 and 7 = 0. Then (JH) and (jl} become 

A' = {py+-y + ^,B = -, 



P 



p- p^' p 
p^ pp 



A^ AA. 



A = 4^ + 2^. (5) 

A^ AA ^ ^ 

Here we explain p as the sum of the density parameters of cold dark matter and baryons, and 
A as the dark energy density. We can designate the corresponding dimensionless densities 
as fip and I^a, respectively, with fip + Da = 1. Then, from we obtain 
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Because the 3D space is flat (A; = 0), we assume that both fi{t) and z/(t) can be expressed as 
power series of t. Using the constraint that A does not decay faster than p does, we consider 
the following choice for fi{t) and z^(t): 



p{t) = t + -,u{t) = r, 



(7) 



where n is a parameter to be determined latter. The corresponding form of A{t, y) is then 
given by 



A' 



1 2 



+ 



1 



(t + a/t) 



2 ' 



tj t + a/t 

If we only consider the condition n > 2 and the late time of the universe t ^ 1, then 
gives 



(8) 



A' = t 



2n-2 



l + O t 



L-2 



j.n—3 



B = {n-l)t 
Substituting ((7j) and Q into (0) and (0), we get 



1 + o {r') 



(9) 



and 



P 
A 



n 



n 



4-2n 



n 



1 + o [r^ 
1 + o (r^ 



(10) 



2{n - 2) 



n + 1 



^ + o{t-')],n..^^^[i^o{t-')]. (11) 

Thus we see that at the late time of the universe t ^ 1, p and A decay with the same 
rate. This is a very useful result in the simulation to the present accelerating universe. 
Remember that we have used the constraint n > 2 in (jU)), (fTUI) and (fTT|. So if n = 2.5, 
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FIG. 1: Evolution of the scale factor ^(t) with = [(t + + + (t + l/t)-2_ There 

is a minimum at t = ife which corresponds to a bounce. 

we have Qp^2/Q ^ 0.22 and Q^jy ^ 7/9 ^ 0.78. If = 3, we have 1/3 ^ 0.33 and 

Q\ 2/3 ~ 0.67. Current observation is about Qp ~ 0.3 and ^ 0.7. So the parameter 
n lies between 2.5 and 3. For an illustration we let n = 3. Then (jH)) gives an exact solution 
being 



A' 



a\ t 



a/t 



+ 



1 



{t + a/t) 



2 ' 



B 



A 



For t 3> 1, this exact solution gives 



(12) 



A 



l + 0{t 

i + oit 



1-2 



, A = 2t~2 
2 



1 + (r^ 
i + ofr^ 



1 + t 



(13) 



Figure 1 is a plot of the scale factor A {t, y) of the exact solution p2|) on the y = 1 hypersur- 
face (with a = 1). From this figure we can see that there is a finite minimum for the scale 
factor A (t) at t = tf, which represents a "big bounce" . Before it the universe contracts, and 
after it the universe expands. 

To see clearly how the universe evolves, we need an appropriate definition for the Hubble 
and the deceleration parameters. Consider the metric (^. On a given y = const hypersurface 
the proper time can be defined as dr = Bdt. So the invariant definitions for the Hubble and 



deceleration parameters should be given as^^ 



1 dA 


1 A 


^^ 


Adr ~ 


BA ~ 


A 


d^A 
^dr^ 


/if 


)■ 







q{t,y) ^ -A^-^/ [^-j^] =-^. (14) 



Then, with use of (jTj) and (fT^ . equation (fT^ yields, for t ^ 1, 



H 



l + 0(t 



1 + t 



L-2 



(15) 



Thus we find that at the late time of the universe t ^ 1, the universe (jl2p is accelerating 
with q ^ —1/2. 

Using (fTB|) in the definition of the proper time dr = Bdt, we find 



T-n^2t (16) 

for t ^ 1 and r ^ r^, where Tf, represents the initial bounce time. So the five- dimensional 
line element gives 

dS^ —.dr'-^ir- n)' (dr^ + r^d^f) - dy^ (17) 
16 ^ ^ 

for r ^ Tf,. This is an approximate metric to describe our present accelerating universe. 

For the special exact solution |T2|l . we can obtain an exact expression for the deceleration 
parameter q{t, y) on a given y = const hypersurface by substituting and (fT^ in the 
second equation of ()14j) . This is a long expression and we plot it in Figure 2 with a = 1 and 
y = l. 

Similarly, we can use the exact solutions p2j) and (|7j) in the definition © to obtain exact 
expressions for Qp and ^a- The global evolutions of Qp and are ploted in Figure 3 with 
a = 1 and y = 1. 

From Fig. 2 and Fig. 3 we can see clearly that the global evolution of the universe is 
divided into four stages separated by the "big bounce" point tb and two critical points td 
and tc2- The first stage is at < t < in which the universe contracts from a A -dominated 
vacuum. The bounce point is a matter singularity (see also Refs. [10,15,21]), across 



7 



4- 


n 








3- 










2- 










1- 




\ t, 













3 


4 5 6 


_-| - 










-2- 










-3- 










-A- 











FIG. 2: Global evolution of the deceleration parameter q(t, y) of the solution (12) with a = 1 and 
y = 1. The bouncing time \s ai t = ti,. There is a critical time td before which the universe is 
decelerating and after which the universe is accelerating. 
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FIG. 3: Global evolutions of J7p and J^a of the solution (12) with a = 1 and y = I. The solid line 
represents Qp and the dashed line represents Oa- There is a critical time tc2 at which 0,\ takes 
over Sip and dominates the universe. 

which flp jumps from — oo to +oo and fl\ jumps from +oo to — oo. The second stage is 
at tb < t < tci in which flp decreases and Oa increases with flp > fl\ and g > 0, and 
the universe is decelerating. The third stage is at td < t < tc2 in which the deceleration 
parameter q becomes negative, implying that the expansion of the universe turns to speed 
up. The fourth stage is at tc2 < t in which < 0, so the universe is accelerating. In this 
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stage, r^A overtakes Qp and the universe is dominated by dark energy. For t ^ 1 we have 
Qa/^p — ^ 2, in agreement approximately with observations of the present stage of our 
universe. 

III. DISCUSSION 

The "big bounce" solution is characterized by having a "bounce", rather than a 
"bang" , as the "beginning" of our expanding universe, and by having a evolving cosmolog- 
ical "constant". Mathematically, the general solution (Q) contains two arbitrary functions 
fi{t) and Different choices of these two functions may give different models to de- 

scribe different stages of our universe. In this paper we find that a simple choice of them 
yields a simple exact solution (fT^ which can describe the present accelerating universe in a 
satisfactory manner. 

We should emphasize that although this simple exact solution exhibits many interesting 
features of the bounce model, it has to be generalized to meet more observations and to 
explain the bouncing from a physical point of view. For instance, the two arbitrary functions 
fi{t) and might be not as simple as given in ((7j) for a real cosmological model. Meanwhile, 
a variable cosmological "constant" might be also too simple to describe dark energy. It is 
known from literatures^'^'^^"^^ that a properly chosen scalar field might be more suitable to 
describe dark energy. Be aware that our model is 5D empty and 4D sourceful with an 
induced 4D energy-momentum tensor which has been supposed to contain a cosmological 
term A as shown in Q. Thus, instead of a A term, we probably should consider the case 
where the induced matter contains a scalar field as a component. We wish this scalar field 
may provide us with a mechanics to explain the acceleration as well as the bouncing of the 
universe. We are going to do this in future studies. 

We thank Baorong Chang, Huanying Liu and Feng Luo for discussion. This work was 
supported by NSF of P. R. China under Grant 10273004. 
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